Bi-partite entanglement entropy in massive two-dimensional quantum field theory 
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Recently, Cardy, Castro Alvaredo and the author obtained the first exponential correction to 
saturation of the bi-partite entanglement entropy at large region length, in massive two-dimensional 
integrable quantum field theory. It only depends on the particle content of the model, and not on 
the way particles scatter. Based on general analyticity arguments for form factors, we propose that 
this result is universal, and holds for any massive two-dimensional model (also out of integrability) . 
We suggest a link of this result with counting pair creations far in the past. 
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Entanglement is a fundamental characteristic of quan- 
tum mechanics: a measurement at a point in space 
may affect instantaneously measurements performed else- 
where in a way that cannot be described by local vari- 
ables. This leads for instance to enhanced performances 
of quantum algorithms. Many measures of quantum 
entanglement between local degrees of freedom in pure 
states are available, but one that has been widely stud- 
ied recently is the entanglement entropy QJ (it appeared 
first in the context of black hole entropy [2]). It measures 
the entanglement between two sectors of compatible ob- 
servables, and is a powerful tool for assessing properties 
of ground states. Its study has already led to a further 
understanding of quantum systems more generally. 

We study the universal behaviour of the entanglement 
between two spacial regions, A and its complement A, 
in the ground state of any extended one-dimensional 
systems. The Hilbert space is H = Ha ® Ha, an d 
the entanglement entropy is the von Neumann entropy 
of the reduced density matrix associated to A: Sa — 
-TT nA p A \og(p A ), PA = Tr w Jgs)(gs| where |gs) is the 
ground state. For a segment A of length R, the universal 
scaling limit, near to a quantum critical point, is obtained 
by making the correlation lengths £ Q , a — 1 , . . . , I and R 
go to infinity with R/£ a = rm a . This defines the dimen- 
sionful length r and a set of particle masses m a (with 
mi < m a \/a) of a quantum field theory (QFT). We will 
assume space-time translation and relativistic invariance, 
with discrete mass spectrum (as occurs in most cases). 

When the length r becomes much larger than all to" 1 , 
the entanglement entropy is known to saturate to a con- 
stant that scales like | log(£i) + U, where c is the central 
charge of the critical point and U is a finite number. Be- 
yond this saturation constant, one expects exponential 
corrections controlled by the masses. However, the exact 
form of these corrections is expected a priori, like the 
constant U and the correlation lengths, to depend in a 
very non-trivial way on the details of the QFT. 

In 0], then in a slightly more general situation in 
the leading correction to saturation of the entanglement 
entropy at large r was evaluated in massive integrable 
QFT, and found to take the same form independently of 



the (factorised) scattering matrix: 



SU = -log(£i) + L/- 



-^^ (2rm Q ) + O(e- 3 ^) (1) 



where Kq is the modified Bessel function. 

In this letter we propose that the same leading cor- 
rection holds generally in massive QFT, not necessarily 
integrable. This is a very general, universal and surpris- 
ing result, as no details of the scattering matrix are in- 
volved, and it only depends on the particle spectrum (a 
basic QFT property). It is also an exact low-energy re- 
sult out of integrability. This is both notoriously difficult 
to obtain (it is impossible to calculate such exponential 
corrections for correlation functions out of integrabiliy) , 
and closer to experimental situations than perturbative 
results. Indeed, because of asymptotic freedom, pertur- 
bation theory of renormalisable models only reliably pro- 
vide high-energy information, whereas most experiments 
on solid state physics occur at energies much below the 
energy gap (this is at the basis of the success of mas- 
sive integrable models, see e.g. [H}). Also, it has been 
observed in integrable models that leading large-distance 
corrections in general give good approximations up to 
very small distances, typically rm\ =0.1. Likewise, for- 
mula ([1]) should give numerically a good approximation 
to the universal entanglement entropy in general QFT. 
At a deeper level, the leading correction contains "clean" 
information about the particle spectrum, an important 
property of a QFT model. It gives all the masses with 
m a < 3to,i/2: "half" of the spectrum, since in non- 
integrable models m a < 2miVa (otherwise the particle 
would decay). 

We provide a derivation of this result based on the scat- 
tering theory of QFT, which we review, and a physical 
interpretation that explains some of its features. 

Hilbert space and local fields. 

The Hilbert space of massive two-dimensional relativis- 
tic QFT is formed by asymptotic states [l^]. There are 
two natural orthonormal bases: particles in the far past 
(in), and in the far future (out). They will be denoted 
by \6 X , ... , Ok)™'™* where 0j, i = 1, . . . , k are rapidi- 



2 



ties and on, i = l,...,k are quantum numbers (with 
vacuum |vac)). In order for particles coming from the 
past to interact, they must have decreasing rapidities 
from left to right. Hence we will adopt the ordering 
01 > . . . > 6k for in states, and the opposite for out 
states. The bases are fixed once impact parameters are 
chosen: extrapolated trajectories meet, say, at x = 0, t = 
0. The energy and momentum are m Qi cosh(0i) and 
X^i m "i sinh(0i), where m a are the masses. 

The Hilbert space alone is not enough to fix a QFT 
model. We need to specify observables. The most im- 
portant is the energy field e(x), with H = j_ dxe(x) 
where x is the position. Local energy measurements 
are quantum mechanically independent at space-like dis- 
tances: [e(x),e(x')] = for x ^ x'. In QFT, we look for 
all local fields (D(x), defined by [0{x), e(x')] = 0, x ^ x' . 

Another observable quantity is the scattering matrix 
(or S- matrix): the amplitude of probability for a given 
in configuration to end up as a given out configuration. 
The LSZ formula gives it in terms of correlation functions 
of local fields ^ a (x) "associated" to the particles. But 
the iS-matrix is believed to completely fix a QFT model, 
hence the opposite also holds: local fields can be deduced 
from it. Below we take this viewpoint. 

The replica trick and branch-point twist fields 

We review the arguments of [3J, |4[. There, the well- 
known replica trick Sa = — lim n _>i -j-Tr^ (Pa) was 
used to evaluate the entanglement entropy. The quantity 
Ti-Ha(Pa) * s identified in the scaling limit as the QFT 
partition function on a Riemann surface with n sheets 
cyclicly connected on the interval A. 

Let us consider instead a model composed of n inde- 
pendent copies of the initial QFT model on R 2 . The 
quantum numbers become doublets (a,j), with j — 
1,2, ...,n representing the copy number, and the lo- 
cal fields acquire an index j. The total energy den- 
sity is e(x) = J2j e j( x )- Particles on different copies do 
not interact (trivial scattering): 101,02)?™ • w ■ \ = 

l^^i)^, i2 ), (ai)jl )forii^i 2 . 

The transformation by which the copy numbers are 
cyclicly permuted is a symmetry of the multi-copy model. 
The branch-point twist field T n is a twist associated to 
this symmetry. It is mainly defined by the equal-time 
exchange relations 

* -(<r'\T (t\ — / Tnix)^ a ,j+l{x ) x' > X , . 

* a > j[X ,Jn[X) -\T n (x)* atj (x') x'<x [Z) 

(with t aj+n = f„j). Other requirements uniquely fix 
this twist field: it is invariant under all other symme- 
tries and has minimal scaling dimension. As n —* 1, it 
becomes the identity operator. 

By symmetry, it is local: [T n (x), e{x')] = (i/ x'). 
Its insertion in a euclidean correlator produces, as func- 
tion of the positions of other fields inserted, a branching 
by which the copies are cyclicly connected through the 



cut extending on its right [l(|. Thanks to this branching, 
the two-point function is proportional to the partition 
function on a Riemann surface, and with an appropriate 
analytic continuation in n, we have 

Sa = flog(£i) + l7- lim A( vac | Tn ( r )T n t(0)|vac) . (3) 
6 n— »i dn 

Below we use locality of branch-point twist fields and 
the iS-matrix of the n-copy model to justify some of their 
properties. 

Form factors of branch-point twist fields. 

Let us consider the simple matrix element 

F^^-e^n) = (vac|T n (O)|0i,0 2 )» !/l2 (0j > 6 2 ) 

(4) 

where /ij are double indices (a*, ji). We used spinlessness 
and relativistic invariance to write it as a function of the 
rapidity difference 9 = 6\ — 62- This gives the function 
for real 9 > only. Let us now consider its analytic 
continuation to complex 9. 

For usual local fields, the analytic structure of such 
matrix elements is well known Q. In terms of Mandel- 
stam's s variable, s — vr? a + m 2 a2 + 2m ai m a2 cosh(0), 
it has a branch cut on s E [(m ai + m Q2 ) 2 , 00), poles for 
s G (0, (m Ql + m a2 ) 2 ) due to bound states, and no other 
singularities on the physical sheet, which covers the com- 
plex plane excluding the cut s £ [(m ai +m Q2 ) 2 , 00). Out 
of integrability, there are other branch points on the cut 
itself due to inelastic scattering, but these do not affect 
the physical sheet. On this sheet, the value on the up- 
per shore of the cut is the form factor with an in state 
as above, and that on the lower shore is a form factor 
with the same particles and momenta, but forming an 
out state. 

It is a simple exercise to translate this analytic struc- 
ture in terms of the variable 9. The strip lm(6) £ (0, 2ir), 
the physical strip, is a double covering of the physical 
sheet (with both shores of the cut separated). The form 
factor with in state is on lm(0) = 0, Ke(9) > and 
lm(6) = 2ir, Re(0) < 0, and that with out state is on 
lm(0) = 0, Re(0) < and lm(0) = 2tt, Re(0) > 0. Poles 
are on lm(0) € (0, 2n), Rc(0) = 0, and are symmetrically 
distributed about the line lm(0) = in. Out of integra- 
bility, = and 9 — 2ni are branch points, since the 
proof that they are ordinary point at integrability relies 
on the lack of particle production. Also there are more 
branch points on the lines lm(0) = and lm(0) = 2n 
(i.e. outside of the physical strip) at inelastic scattering 
thresholds, and in general the analytic structure outside 
of the physical strip on the 0-plane is very complicated. 

An important aspect of the present paper compared 
to previous works at integrability [1, 0| is to show that 
we only need the analytic structure on the physical strip: 
in particular the inelastic-scattering branch points on its 
boundary do not affect the derivation. 

For the branch-point twist fields T n , these properties 
are modified. In order to understand them, we will use 
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FIG. 1: Rotating the 81 wave packet in the euclidean plane. 
A shift by more than 3in/2 brings it from copy j\ to copy 
ji + 1 (we chose the branch cut to go straight to the right). 

an intuitive picture in the euclidean theory, where an 
imaginary shift of rapidity corresponds to a rotation of 
the wave packet (see fig. 1). A complete explanation 
would require the use of Feynman diagrams, but devel- 
oping these details is out of the scope of the present pa- 
per. However the arguments we present show that the 
well-known analytic structure described above for ordi- 
nary fields is a direct consequence of this intuitive picture. 
Also, more general analytic properties of twist fields form 
factors in integrable models were derived following simi- 
lar ideas in combination with integrability, and could be 
verified to a large extent thanks to exact solutions. This 
gives strong support to arguments from this intuitive pic- 
ture in general. Although twist field form factors were 
only studied at integrability until now, the arguments, 
as is clear below, hold also out of integrability for two- 
particle form factors. 

First, the value on 9 < is as usual: 

F WlMa (0i - 9 2 ,n) = (vs,c\T n (0)\e 1 ,6 2 )™% 2 (9 2 > 9,) . 

Intuitively, from particles in an in-state, continuously 
changing 9\ — 9 2 to make it negative gives particles trav- 
elling away from each other, which is an out state. In the 
case where j\ ^ j 2 , particles do not interact, so that 



F, 



(5) 



Second, the value at 9 + 2iri for real 9 can be obtained 
by rotating the wave packet of particle 1 clockwise. As 
the rapidity 9\ arrives at 2iri, particle 1 is on the right 
of particle 2 and they are travelling in the out state con- 
figuration. However, particle 1 is now on copy j\ + 1: it 
is the state \02,Oi)°^n t with fii — + 1 mod n). 

Hence we have 



F„ 



! (9 + 2iri,n) = F fl2jll (-9,n) 



(0) 



Third, the function F^ 1 ^ 2 (9,n) has poles at purely 
imaginary values of 9 for which bound states of particles 
ai and a 2 provide additional on-shell channels. These 
bound states can occur only if the particles travel on the 
same copy when they interact. Shifting Im(6?i) by an 
amount smaller than 7r, particles 1 and 2 interact on the 
left of the field T n , where they are respectively on copies 
ji and j 2 . Hence for lm(0) <E (0, w), poles may occur only 
for ji = j 2 . On the other hand, going beyond n, they 
interact on the right, where they are on copies j% and 
j 2 — 1, or ji + 1 and j 2 . Hence for lm(9) S (n, 2%), poles 
occur only for j\ = j 2 — 1 mod n. 



Finally, there is an additional pole at 9 = itr, a "kine- 
matic pole" . Consider d\ — > 6% + in. The resulting wave 
packet represents a single-particle out state with (real) 
rapidity 9i and particle type (a\,ji), that is ™ t j 1 (9i\, 
where a± is the anti-particle of ct\ . In regions of space 
where the in and out wave packets correspond to the 
same particle type and the same copy, they overlap. 
There are two distinct regions where overlaps may occur: 
x — > 00 and x — > —00. They provide the main contribu- 
tion to the associated matrix element for 9\ ~ 9 2 : 



fei ,0, {T n )m Vcosh 0! cosh 9 2 J dxe -«™(^he 1 -sinhe 2 ) x 
x (^i,i20(-aO + s n+h32 e(x)) 

where Q(x) is the step function, and Kronecker delta 
functions of copy numbers are modulo n. This gives both 
a Dirac delta function at 9\ =62, and a pole (its principal 
value) as 9\ — > 9 2 . Only the pole can be seen from the 
analytic continuation in 9\\ 



iF, 



1 + in) 



(Tn) 



( S h,j2 - 5 h+l,j 2 ) ■ ( 7 ) 



Apart from all these poles, the form factors are analytic 
in the physical strip lm(9) € (0, 27r). Using |5j) and (J6j) , 

f («ui),(a2j 2 )( 9 ) = F {aulh{a2A) (9 + 2ni(j 1 - j 2 )) (8) 

for < ji — j% < n — 1. Then, from the pole condi- 
tions and ([7]), -F( Qi ,i).(q 2 ,i)(#) is analytic in the extended 
physical strip In\{9) € (0, 27m) except for possible bound- 
state poles at purely imaginary 9 on Im(#) S (0, 7t) and 
lm(0) G (27m — 7T, 27m), and for kinematic poles, when 
ai = a 2 , at 9 = in and 9 = (2n— l)i7r with residues i(T n ) 
and —i(T n ). This is the structure found in 0,0] in inte- 
grable QFT. The larger analytic region of i 7, ( Qll i).(Q 2 ,i)(^) 
is surprising, since for usual local fields in non-integrable 
QFT there are branch points on Im(#) = 2ir. It was 
instrumental that particles on different copies do not in- 
teract, and that a shift by 2?ti changes the copy number. 
Beyond the extended physical strip, the analytic struc- 
ture is much more complicated. 
Entanglement entropy. 

The two-point function in ([3]) can be expanded at large 
distances by inserting intermediate states between the 
twist fields and using space translation covariance. In 
particular, the two-particle contribution is: 



x e 



d8id9 2 f a ^(9 1 -8 2 ,n) x 

r(m a cosh 0\ cosh O2 ) 



(T n ) 2 f a ,f3(9, n) = J2 |*W)(/3,i)(0 + 2^j,n)| 2 (10) 

3=0 

In order to evaluate the entanglement entropy, we need to 
analytically continue from n = 1,2,3,... to n € [1, oo), 
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then to take the derivative at n = 1. This was done in 
[H, 0] in various ways. We recall here the main steps, em- 
phasizing that we only need the analytic structure on the 
extended physical sheet, independently from the struc- 
ture on its boundary (where there are inelastic-threshold 
branch points) or further away. 

First for n £ [1, oo), Ft a A)(p,i){9) has poles for a = j3 
at 9 — in and 9 = 2iim — m with unchanged residues, 
and no poles in the strip lm(9) G (n, 2nn — ir). As noted 
in 0], this follows from the picture of a space with two 
conical singularities of angle 2nn, valid for real positive 
n: kinematic poles for particles going past the conical 
singularity on its left and right, and no bound states in 
the extra space. Second, form factors vanish as n — > 1 
since %,{x) — > 1. A vanishing like n — 1 was observed 
in 0, |j| j an d we assume this still holds. This implies 
no contribution from one-particle form factors. The two- 
particle contribution comes from the non-uniform con- 
vergence of form factors as n — > 1, due to the colli- 
sion of kinematic poles, as explained in More pre- 
cisley, the sum in (|T0l) can be done by Y^j=a f(j) = 
/(0) + (2*)- 1 § dzf(z) cot nz - Q = /(0) + f™ x dy(f(n + 
iy — 7) cot n(iy — 7) — f(iy-\-j) cot n(iy + 7))/2 — Q with 
< 7 < 1/2 and where Q cancels the residues of f(z). 
Note that the boundary of the extended physical sheet 
is avoided. Then, as n — > 1 the kinematic residues in 

Q collide and give [jj (^/o,/j(^»)) n=; = T , *W*a lj 9 (a 
formula like this was obtained earlier in the free Dirac 
model in 0). This gives (TJ). 
Discussion: entangled pairs. 

The most striking feature of fl} is its independence 
from the scattering matrix. If the entanglement entropy 
Sa were counting a "number of links" connecting points 
of A and A, we would have (a regularised version of) 
Sa = f A dx Ja dx's(x — x') for any region A. This agrees 
with Sa = Sz and holds explicitly in some valence-bond 
descriptions is equivalent to 

saying that the mutual information of [lj| is extensive, 
or to Sa = - J2 X <x'edA r l{ x ) r l(. x ')S{x,x'\ where rj(x) is ±1 
if x is a left /right boundary of a connected component. 
However, there is convincing numerical evidence that the 
latter formula fails (although slightly) in general two- 
dimensional QFT [l3|. Yet, factorisation of twist-field 
correlation functions immeditaly implies extensivity of 
the mutual information in the large-distance limit. Hence 
the failing should only be due to some "non-locality" 
at the end-points of the "links" , and we may interpret 

s(x) = — \ d S ^°i^ 1 for x much larger than the correlation 
length, as measuring an "entanglement density" , a corre- 
lation between quantum disturances a distance x apart. 

A natural candidate is from virtual pairs created far in 
the past (see fig. 2), in a similar way to what happens 
in quenching [l4J . Technically, the collision of kinematic 
poles indeed says that we are considering particles com- 
ing from a common virtual pair. The probability that an 



x' x 




FIG. 2: A pair contributing to the entanglement density. 

entangled pair survives for a time t is essentially ruled by 
quantum uncertainty principles, cx e~ tE where E is the 
total energy, independently from the scattering matrix. 
The trajectories are linear on the world sheet since in one 
dimension, conservation of momentum and energy for- 
bids smooth changes of their rapidities. Hence the main 
contribution to s(x) at large distances should be indepen- 
dent from the scattering matrix. For instance, a pair of 
particles of mass m and velocities v and —v has life-time 
E^ 1 = (2m,) -1 vl — v 2 . When hitting points x and x' , 
it survived for a time t — (x — x')/(2v). Its contribution 
is m 2 dv f(v)e~ tE , which should be positive (f(v) > 0), 

and s(x — x') ~ J2a=i m a Jo ^ v f( v ) e " %/l ~" 2 ■ The ex- 
ponent is maximum at v = 1/V2, where it has value 
—2m a (x — x'), in agreement with the correction terms 
(HJ). Also, UJ implies f(v) = 32v ^l_ v '2yi , which indeed is 
positive for v € (0, 1). 

In conclusion, we evaluated the leading large-distance 
correction terms for the entanglement entropy in any uni- 
tary two-dimensional QFT. This is one of the few ex- 
amples of an exact low-energy result, and of the use of 
the analytic structure of form factors, out of integrabil- 
ity. The heuristic arguments we provided give a physi- 
cal explanation of the main features: independence from 
the scattering matrix, and why the collision of kinematic 
poles lead to this result. It would be useful to have a 
derivation of the form factor properties from perturba- 
tion theory, a better study of the analytic continuation 
in n involved, and a better understanding of extensivity 
properties of the mutual information. 
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